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Abstract 


In recent years, fault tolerant flight control systems have gained an increased 
interest for high performance military aircraft as well as civil aircraft. Fault 
tolerant control systems can be described as either active or passive. An active 
fault tolerant control system has to either reconfigure or adapt the controller in 
response to a failure. One approach is to reconfigure the controller based upon 
detection and identification of the failure. Another approach is to use direct 
adaptive control to adjust the controller without explicitly identifying the fail- 
ure. In contrast, a passive fault tolerant control system uses a fixed controller 
which achieves acceptable performance for a presumed set of failures. 

We have obtained a passive fault tolerant flight control law for the F/A-18 
aircraft which achieves acceptable handling qualities for a class of control sur- 
face failures. The class of failures includes the symmetric failure of any one 
control surface being stuck at its trim value. A comparison was made of an 
eigenstructure assignment gain designed for the unfailed aircraft with a fault 
tolerant multiobjective optimization gain. We have shown that time responses 
for the unfailed aircraft using the eigenstructure assignment gain and the fault 
tolerant gain are identical. Furthermore, the fault tolerant gain achieves MIL- 
F-8785C specifications for all failure conditions. 
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1 Introduction 


1.1 Background 

Aircraft flight control systems are designed with extensive redundancy to ensure a 
low probability of failure. During recent years, however, several aircraft have expe- 
rienced major control system failures. These have caused an increased interest in 
fault tolerant flight control systems. The objective of fault tolerant flight control 
is to control and safely land the aircraft in case of severely damaged or inoperable 
control surfaces. The two approaches to fault tolerant control are active and passive 
control. An active fault tolerant control system has to either reconfigure or adapt 
the controller in response to the failure. One way is to reconfigure the controller 
based upon detection and identification of the failure. Another way is to use direct 
adaptive control to adjust the controller without explicitly identifying the failure. In 
contrast, a passive fault tolerant control system uses a fixed controller which achieves 
acceptable handling qualities for a given set of failures. 

Several authors have utilized eigenstructure assignment to design reconfigurable flight 
control systems. Gavito and Collins [1] used eigenstructure assignment to recover the 
undamaged modal response under the assumption that the failure has been detected 
and identified. Napolitano and Swaim [2] used eigenstructure assignment to remove 
the lateral-longitudinal coupling induced by an asymmetric control surface failure. 
Jiang [3] used eigenstructure assignment to recover the dominant eigenvalues and 
eigenvectors of an aircraft longitudinal control system which has undergone some 
operating condition variations or system component failures. However, all of these 
approaches require an on-line identification of the parameters of the plant after failure. 

Jiang and Zhao [4] used eigenstructure assignment to design passive fault tolerant 
controllers for actuator failures. This approach is based on robust regional eigenvalue 
assignment with the addition of a precompensator. However, this design method re- 
quires extensive redundancy. Gorrec et. al. [5] proposed a multimodel approach for 
the landing phase of a large transport aircraft. This method uses a bank of models 
covering the entire flight envelope. First, the authors design an initial gain for a cho- 
sen original model. This gain is used to control all of the models and a multimodel 
analysis is used to detect the worst model. Then, a quadratic optimization procedure 
is used to improve the behavior of the worst model while keeping good performance 
relative to the original model. However, the designer must be careful to avoid con- 
flicting objectives. 

We have obtained a new result for eigenstructure invariance when it is known which 
surface is most likely to fail. We seek a constant gain output feedback controller such 
that the dominant eigenstructure is invariant under this failure. We show mathemat- 
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ically that the solution is that the failed surface should not be used. That is, if the 
j th control surface will fail, then the j th row of the constant output feedback gain 
matrix will be simply zero. For this purpose, first we derive a basis for the subspace 
in which the eigenvectors of the failed system must lie. Then, we show that there 
exists a constant output feedback gain matrix which yields an invariant dominant 
eigenstructure both before and after failure. Finally, we prove that the j th row of this 
feedback gain matrix will be zero. 

We have obtained a passive fault tolerant flight control for a control surface failure 
when it is not known in advance which surface will fail. We consider the linearized 
lateral dynamics of F/A-18 aircraft. The class of control surface failures we consider 
includes the symmetric failure of any one control surface being stuck at its trim value. 
It is not known in advance which surface will fail. We computed an optimal feedback 
gain using an off-line multiobjective optimization technique. A comparison was made 
of an eigenstructure assignment gain designed for the unfailed aircraft with the fault 
tolerant multiobjective optimization gain. We have shown that the time responses 
for the unfailed aircraft using the eigenstructure assignment gain and the fault tol- 
erant gain are identical. Furthermore, the fault tolerant gain achieves MIL-F-8785C 
specifications for all failure conditions. 
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1.2 Contributions of the Work 

• We have a new result for eigenstructure invariance when it is known which surface 
is most likely to fail. We derived a basis for the subspace in which the eigenvectors of 
the failed system must lie. We showed mathematically that, if the j th control surface 
will fail, then the j th row of the constant output feedback gain matrix will be simply 

zero. 

• We have obtained a passive fault tolerant flight control for a control surface failure 
when it is not known in advance which surface will fail. We used an off-line multi- 
objective optimization technique to design an optimal control for a predefined class 
of control surface failures. We have shown that the time responses for the unfailed 
aircraft using the eigenstructure assignment gain and the fault tolerant gain are iden- 
tical. Furthermore, the fault tolerant gain achieves MIL-F-8785C specifcations for all 
failure conditions. 


2 Application of Eigenstructure Assignment to Flight 
Control Design 

2.1 Literature Review 

Consider an aircraft modelled by the linear time invariant matrix differential equation 
described by 

x = Ax 4- Bu 0) 

y = Cx ( 2 ) 

where x is the state vector ( n x 1), ti the control vector (m x 1) and y the output 
vector (r x 1). It is assumed that the m inputs and the r outputs are independent. 

If there are no pilot commands, the feedback control vector u equals a matrix times 
the ouput vector y: 

u = -Fy ( 3 ) 

Theorem [6]: 

Given the controllable and observable system described by Eqs. 1 and 2. max(m.r) 
closed loop eigenvalues can be assigned and max{m,r ) eigenvectors (or reciprocal 
vectors by duality) can be assigned with min(m,r) entries in each vector arbitrarily 
chosen using constant output feedback. 

Andry et. al. [7] concluded that the eigenvectors v t must be in the subspace spanned 
by the columns of (A, 7 - A)~ } B. This subspace is of dimension m which is equal to 
the number of independent control variables. Thus, if we choose an eigenvector r, 
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which lies precisely in the subspace spanned by the columns of (A J — A) l B, it will 
be achieved exactly. In general, however, a desired eigenvector vf will not reside in 
the prescribed subspace and hence cannot be achieved. 


Ref [7] gives a way to find the ” best possible choice ” for an achievable eigenvector. 
This best possible eigenvector is the projection of vf onto the subspace spanned by 
the columns of (A,/- A)~ l B (in the least square sense). In many practical situations, 
complete specification of v t is neither required nor known, but rather the designer is 
interested onlv in certain elements of the eigenvector. Thus, assume that V{ has the 
following structure: 

vf = [Uii,X,X,X,X,Vy,X,X,Vi„] T 

where v i3 are designer specified components and x is an unspecified component. De- 
fine, as shown by Andrv et. al. [7], a reordering operation { }*' such that 


{vf} R ' = 


Li 

d{ 


(4) 


where U is a vector of specified components of vf and d* is a vector of unspecified 
components of vf. The rows of the matrix (A J - A)~ l B are also reordered to conform 
with the reordered components of vf. Thus, 


{(A,/ - A)~ l B} Rt = 


U 

A 


(5) 


Then, as shown by Andry et. al. [7], the acheivable eigenvector vf is given by 
vf = (X i I-A)- 1 Bz i (6) 

where Zi = L\li and where (-) 1 denotes the appropriate pseudoinverse of (•). 


2.2 Problem Statement 

We consider the linearized lateral dynamics of the F/A-18A aircraft. The rigid body 
states are lateral(side) velocity(z;), yaw rate (r), roll rate (p), and bank angle (0). 
These states are augmented with first order actuator dynamics and a yaw rate washout 
filter to yield a 9 th order model. The control surfaces are asymmetric stabilator (<5 SC ), 
asymmetric trailing edge flaps ( <5( ec ), ailerons (^ac)i an d rudder (<5 rc ). The inputs are 
stabilator command (c^ c ), trailing edge flaps command (6 tec ), ailerons command (<5 ac ), 
and rudder command (<5 rc ). The measurements are sideslip angle (/?), washed out yaw 
rate (r„, 0 ), and roll rate (p). A control surface failure is modeled by setting the cor- 
responding column in the control derivative matrix to zero. Each failure condition 
corresponds to one control surface being fixed at its trim value. 
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The aircraft can be described by 


x = Ax + Bu 

y - Cx 

where dim [x\ = n , dim [u] = m, dim [y] = r, and where 


( 7 ) 

( 8 ) 


A = 


B t = 
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Observe that Eqs. (9)-(ll) may be rewritten compactly as 


A = 


B 

C 


= 6 


dim 0 mx (n-m) 

A21 A 22 

In, 

0 (n— m)xm J 


= 0. 


(C 12) rx(n-m) J 


( 12 ) 

(13) 

(14) 


and the aircraft with a control surface failure may be described by 

x = AfX + Bu 

y = Cx 


(15) 

(16) 



where 


( 17 ) 

( 18 ) 


A, 

A.4 

Ay^i 


A + AA 

0 m 0?m x (n-m) 

A 7^21 0( n -- rn )x(n-m) 

0 444 &(m+l)ji * 0 

0 * 0 


( 19 ) 


2.3 Eigenstructure Invariance for an Anticipated Flight Con- 
trol Surface Failure 

Consider an LTI plant which is augmented with first order actuator dynamics. Sup- 
pose we know which control surface is most likely to fail. We seek a constant gain 
output feedback controller such that the dominant eigenstructure is invariant under 
this failure. We will show mathematically that the solution is that the surface which 
will fail should not be used. That is, if the j th control surface will fail, then the 
j th row of the constant output feedback gain matrix will be zero. We present some 
interesting intermediate results. First we derive a basis for the subspace in which 
the eigenvectors of the failed system must lie. Then we show that a constant output 
feedback gain matrix exists which yields an invariant dominant eigenstructure both 
before and after failure. Finally we show that the j th row of this feedback gain matrix 
will be zero. 
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Lemma 1 : 

Let Af, i = 1, 2, • ■ • , r be the set of desired eigenvalues. 

A basis for the m-dimensional subspace spanned by the columns of (Af / - .4 
given by 



' 0 ' 


a\S 


0 


' 0 ' 

0 






GjS 




, , . 

J 

0 


0 

0 




Oim^ 

. 7i . 


0 


. 7m . 


. 7; . 





where -a* = 

and 7 = J 7i> ' • • > 7m ] = &{Adj{A 1 )) 21 

Proof 01 Lemma 1: 


(Afz-Ay) = 


(Af + S)I m Omx(n-m) 
—^21 + A^i Af/ — A22 


V41-A,)-' 


Adj(\fl - Af) 
det(\fl — A/) 


Let 

(Af/ - Aj) = A 1 


Let Adj{A l ) = [c pq } T where the cofactors c vq are given by 
c vq = (-1 Y + « ■ det{Al q ) 

Due to the special structure of A 1 , it follows that 


(Af -I- S) ■ det(A l pq ) if p = q, P, q = 1, 2, • • 
0 ifp^g, p,q= 1,2, •• 


m 

m 


Therefore, the adjoint of A 1 may be written as 


Adj(A') = 


(Af + <5)r (Ad^A 1 ))^ 

(Adj(A 1 )) 21 (Adj(A ] )) 22 


;)-'// is 

( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 

(26) 
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Where T = diag{det(A l n ), det(A\ 2 ), • • • , detiA^J) 
Finally 


(A «I-A f )- l B 


Adj(A^) • B _ 1 _ 

det(A l ) det(A 1 ) 


(xf + <5)<5r 

6(Adj(A l )) 2 i 


Define 


OL — 


( Af + .5)<sr 

det{A l ) 


(\*+6)6det(A{ lc ). n 
— HeHAU 1 'A — A? ! 


m 


a * — dei(A l , 

7=[ ] =^(A! , ))2, 

Then, the basis vectors are given by (20). 


Theorem 1 (Existence): 

There exists a constant output feedback controller u = Fy such that: 
(Af + BFC)v( = (A + BFC)v( = Aft// ; i = 1, 2, • • • , r 

Proof of Theorem 1 : 

Use eq.(15) to obtain 

(j4/ + BFC)vf = (A + AA + BFC)v{ = (A + BFC)v{ + AAv/ 
Let the j th control surface fail. Then 


AAv{ = 



0 TO 

a (m+l)j 


0 ’ Clnj 


1 

? 

1 

£ 

X 

£ 

o 

o 


1 

^ .. 
1 

: 0(n— m)x(n— m) 

0 


1 

sV' 

1 


From lemma 1, v{ may be written as 


v[ - (3i 


ot\6 

0 


7i J 


+ ■■• + {3 j 


0 

0 


Aj 


+ Pr> 


0 

QZmd 


Am J 


(27) 

(28) 
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(30) 

(31) 

(32) 


(33) 
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Observe that all the basis vectors have a zero in row j with the exception of basis 
vector j . 

Choose fy = 0. Then, AAv{ = 0 and (30) is proven. 

Theorem 2(Gain Computation): 

The constant feedback gain F from Theorem 1 has all zeros in row j. 

Proof of Theorem 2 : 

Rewriting the eigenvalue-eigenvector equation we obtain 

(X t I - A f ) v{ = BFCvf (34) 

We partition (34) conformally mindful of the special structure of Af,B matrices. 


(A* + 8)I m 

0mx(n— m) 


' Z, ' 

= s 

Im 

FC 

' Z, ‘ 

— ^21 + A^21 

Aj/( n -ro) — ^22 


Wi 


0 (n-m)xm 


H’,- 


Where 


v{ = 


Zi 


(35) 


(36) 


[ (A, + S)I m Omx(n-m) 

(At + = 8FCv { ; i = 1 ? • • * , t 



z 


I 


r Zi 1 


z. 


Wi 

= 6 

1 m 

FC 

w t . 

= SFC 

h; 

- 


Omx (n— m) 


(37) 
(38) 


Thus 


/ Aj + 6 

1 5 


)Z t 


Let 


FCv { ; i = 1, 


Z = 



^)Zt, (^)Zr] 

v' = [ v{ , v{, ■ • ■ , v{ j 


Therefore 


Z = FCv } 

and 

F = Z{Cv f )~ l 


(39) 

(40) 

(41) 

(42) 

(43) 
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Finally, 


x I 

where the zero is in row j; j < m (44) 

L x J 

So 

F= 0 0 ••• 0 (CVf)' 1 (45) 

where the zeros are in row j of matrix Z. Then the result follows. 
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Example 

We use the plant described in Section 2.2. We consider the rudder failure which is 
modeled by setting the 4 th column of the control derivative matrix to zero. 

First, we design an eigenstructure assignment controller for the unfailed condition by 
assigning the dutch roll mode to have a damping of 0.707 and a natural frequency of 
2.83 r/s. The roll subsidence mode is assigned to its open loop value of -2. 76. The 
desired dutch roll eigenvectors are chosen for a sideslip and yaw rate mode which is 
decoupled from the roll rate and bank angle. The desired roll subsidence eigenvector 
is chosen for a roll rate mode which is decoupled from vertical velocity and yaw rate. 
Table 1 shows the output feedback matrix computed for the unfailed aircraft . The 
initial condition responses using the eigenstructure assignment gain computed for the 
no failure condition are shown in Figure 1. Observe the degradation in the response 
after failure. 

Next, we compute an optimal feedback gain using Theorem 2 by imposing simulta- 
neous specifications on the unfailed and rudder failure conditions. We require the 
dutch roll and roll modes for both the unfailed and failed conditions to be the same 
as for the eigenstructure assignment controller. Furthermore, we require that the 
fourth row of the dutch roll and roll eigenvectors be zero for both the unfailed and 
failed conditions. Table 2 shows the output feedback gain matrix. Observe that the 
fourth row of this gain matrix is exactly zero. The initial condition responses using 
the optimal gain are shown in Figure 2. Observe that the eigenstructure invariance 
is achieved but with increased coupling to roll rate and bank angle. To reduce the 
coupling, the number of eigenvector objectives is reduced by only requiring the fourth 
row of the dutch roll eigenvector to be zero. The feedback gain matrix is shown in 
Table 3 where we observe that the gains in the fourth row are small but nonzero. The 
initial condition responses using this optimal gain are shown in Figure 3. We observe 
that the coupling to roll rate and bank angle has been significantly reduced, albeit at 
the expense of exact eigenstructure invariance. 


Table 1: Gain matrix computed for no failure 


0 

y\vo 

P 


0.2412 

0.0562 

-0.0009 

stabilator 

0.1317 

0.1742 

0.0038 

trailing edge flap 

0.1352 

0.1918 

0.0047 

aileron 

-0.8366 

0.9019 

0.0379 

rudder 


13 



Table 2: Gain using theorem (dutch roll and roll modes) 


$ Two P 

-9.1724 -5.6198 -0.2616 

4.7829 4.8542 0.1423 

-13.0675 5.2913 0.2041 

0 0 0 


stabilator 
trailing edge flap 
aileron 
rudder 


Table 3: Gain using theorem (dutch roll only) 


0 

^ WO 

P 


8.2678 

-7.6184 

-0.1597 

stabilator 

-4.0159 

2.7116 

0.0564 

trailing edge flap 

-3.9988 

4.5213 

0.0944 

aileron 

0.0613 

0.0298 

0.0198 

rudder 



Gain computed for no failure 




time (sec) 


Figure 1: Gain computed for no failure 
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bank angle (deg) roll rate (deg/sec) ygw rgtg (deg/sec) si ( desli P an 9 le < de 9) 



Figure 2: Gain computed using theorem for dutch roll and roll 
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sideslip angle (deg) 



Figure 3: Gain computed using theorem for dutch roll mode 
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2.4 Passive Fault Tolerant Control for an Unknown Control 
Surface Failure 

We have designed an optimal gain which is tolerant to a control surface failure. The 
class of control surface failures we consider includes the symmetric failure of any one 
control surface with the surface being stuck at its trim value. It is not known in ad- 
vance which surface will fail. We compare an eigenstructure assignment gain designed 
for the unfailed aircraft with a multiobjective optimization gain designed to achieve 
(1) the same response for the unfailed conditon as was obtained using eigenstructure 
assignment and (2) MIL-F-8785C specifications for all failure conditions belonging to 
our class. 

We use program ATTGOAL from the MATLAB Optimization Toolbox [8] to solve the 

multiobjective goal attainment problem described by minimize 7 such that 

x, 7 

F(x) -wj< goal where F(x ), w , and goal are given. 

The function F(x ) is the objective function to be minimized at the point x. ATTGOAL 
attempts to minimize the function values to attain the goal values given by goal. 
Alternatively, to make an objective function as near as possible to a goal value, op- 
tions(15) is set equal to the number of objectives required to be in the neighborhood 
of goal values. These objectives must be partitioned into the first elements of the 
function F(x). Goal is a vector of values that the objectives attempt to attain. The 
weighting vector w is used to control the relative under-attainment or over-attainment 
of the objectives. The weighting function may be set to w = abs(goal) to obtain the 
same percentage of under or over-attainment. When w is positive, ATTGOAL at- 
tempts to make the objectives less than the goal values. To make the objective 
functions greater than the goal values, w is set to be negative. 

We consider the linearized lateral dynamics of the F/A-18A aircraft which was de- 
scribed in Section 2.2. We design an eigenstructure assignment controller for the 
unfailed condition by assigning the dutch roll mode to have a damping of 0.707 and 
a natural frequency of 2.83 r/s. The roll subsidence mode is assigned to its open loop 
value of -2.76. The desired dutch roll eigenvectors are chosen for a sideslip and yaw 
rate mode which is decoupled from the roll rate and bank angle. The desired roll 
subsidence eigenvector is chosen for a roll rate mode which is decoupled from vertical 
velocity and yaw rate. 

Table 4 shows the dutch roll and roll subsidence mode desired eigenvectors, the eigen- 
vectors from the projection onto the achievable subspaces, and the eigenvectors of 
(A + BFC ). First, we remark that the entries in the imaginary part of the dutch roll 
eigenvectors, which we desire to be zero, are not zero after projecting onto the achiev- 
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Table 4: Desired and Achievable Eigenvectors (actuator entries not shown) 


Desired Projection (A + BFC ) 


D. R. 

R. 

D. R. 

R. 

D. R. 

R. 


1 + jx 

0 

1+jO 

0.0001 

1 — j59.4 

0.0001 

r 

xTj 1 

0 

0.0028 + jl 

-0.0168 

0.0032 + jl 

-0.0168 

r 

o + jo 

1 

—0.0137 + j3.83 

1 

—0.012 + j5.24 

1 

p 

0 + jO 

X 

0.0006 + j2.14 

-0.362 

-0.0003 - j2.52 

-0.362 

<p 

x + jx 

0 

-0.0008 + j0.034 

0.0037 

-0.0008 + j0.08 

0.0037 

•r 9 


able subspace. This will cause coupling to the roll rate response. The reason that we 
do not achieve the zero eigenvector entries is that the L { are poorly conditioned. So 
to achieve a gain matrix with reasonable magnitude, we discard the smallest singular 
value when computing the pseudo inverse to obtain z*. If less coupling is required 
then another choice for the desired dutch roll eigenvectors should be made which 
specifies fewer entries. Second, we remark that we have r < m which is different from 
most aerospace applications of eigenstructure assignment which have appeared in the 
literature. Some authors have suggested assigning the left eigenvectors when r < m. 
However, the right eigenvectors still lie in the subspaces spanned by the columns of 
(XJ - A)~ X B. Since we only assign r eigenvalues, we can also assign m entries in the 
corresponding eigenvectors. This is independent of the relation between r and m. 

The feedback gain matrix is shown in Table 5. The sideslip angle responses to a 
one degree initial sideslip for the unfailed and four failure conditions are shown in 
Figure 4 where the oscillitory response corresponds to a rudder failure. We remark 
that the response is unacceptable when a rudder failure occurs. In fact, the dutch 
roll damping is only (, = 0.18 when a rudder failure occurs. This is signficantlv less 
than the minimum value of £ = 0.4 which is required in MIL-F-8785C. We further 
remark that the asymmetric stabilator, asymmetric trailing edge flaps, and ailerons 
are all good producers of rolling moment, but only rudder is a good producer of 
yawing moment. Thus, a rudder failure is the most difficult condition for achieving 
desirable handling qualities. Next, we design an optimized controller by imposing 


Table 5: Eigenstructure Assignment. Feedback Gain Matrix 


P 

T wo 

P 


0.2412 

0.0562 

-0.0009 

stabilator 

0.1317 

0.1742 

0.0038 

trailing edge flap 

0.1352 

0.1918 

0.0047 

aileron 

-0.8366 

0.9019 

0.0379 

rudder 
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simultaneous specifications on the unfailed and four failure condtions. We use the 
multiobjective optimization method which is implemented in the MATLAB Opti- 
mization Toolbox [8] program ATTGOAL. Our objective is for the optimal unfailed 
responses to be the same as the eigenstructure assignment unfailed responses while 
also achieving MIL-F-8785C specifications for the failed responses. The MIL-F-8785C 
specifications for Level 1, Category A (CO and GA), Class IV flight are as follows: 
dutch roll damping and natural frequency greater than 0.4 and 1.0, respectively; roll 
subsidence time constant less than 1 second; and spiral time to double amplitude 
greater than 12 seconds. We remark that the eigenstructure assignment controller 
achieves all these specifications with the exception of the dutch roll damping which 
is 0.18 when a rudder failure occurs. 
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sideslip angle (deg) sideslip angle (deg) 


eigenstructure assignment gain 



Figure 4: Sideslip angle responses 
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We compute an optimal feedback gain using the following objectives: (1) the 
dutch roll and roll subsidence modes are to be the same as for the eigenstructure 
assignment controller, and (2) the dutch roll and roll subsidence modes are to achieve 
MIL-F-8785C specifications for the four failure conditions. We remark that we cannot 
directly specify the spiral mode when using a stability augmentation system (SAS) 
which does not feedback bank angle. 

The optimization parameters are (1) the real part of the dutch roll eigenvalue Re(Ad r )> 
(2) the imaginary part of the dutch roll eigenvalue Im(Ad r ), (3) the roll subsidence 
eigenvalue A ro n, (4) the free eigenvector parameters z ro u for the roll subsidence mode, 
and (5) the free eigenvector parameters z dr for the dutch roll mode where real arith- 
metic is used. The parameters are initialized at the values obtained when using the 
eigenstructure assignment gain. 

The objectives and weightings at the unfailed condition are 

A roll = -2.76 , w = abs(-2.76) 

Re(A dr ) = -2 , w = abs(— 2) 

Im(Adr) = 2 , w = abs(2) 

The objectives and weightings at the 4 failure conditions are 

a >0.4 , w l = -h t = l,...,4 
(w n )d r > 10 , w* = -1; i = 1, ... ,4 

The optimal gain is shown in Table 6 which required 2.04 seconds of CPU time and 
52 function evaluations. The sideslip angle response to a one degree initial sideslip is 
shown in Figure 4. We observe that all responses achieve the desired specifications. 
The dutch roll damping ( dr is now 0.4 for the rudder failure conditon which achieves 
the MIL-F-8785C specification. The yaw rate responses for an initial sideslip for the 


Table 6: Optimal Feedback Gain Matrix 


0 

^WO 

P 


5.2056 

-2.851 

-0.0605 

stabilator 

-2.1993 

1.4108 

0.0292 

trailing edge flap 

-2.651 

1.9134 

0.0400 

aileron 

-0.0766 

0.4808 

0.0293 

rudder 


eigenstructure assignment and multiobjective optimaization gains are shown in Figure 
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5. Again we observe that the eigenstructure assignment gain exhibits an unaccept- 
able oscillation when rudder failure occurs whereas the multiobjective optimization 
gain yields an acceptable response. The roll rate responses for the eigenstructure 
assignment and multiobjective optimization gains are shown in Figure G. Here we 
observe that the roll rate response with the multiobjective optimization gain exhibits 
a significant coupling when a stabilator failure occurs. The bank angle responses 
for the eigenstructure assignment and multiobjective optimization gains are shown in 
Figure 7. Again, we observe that the response with the multiobjective optimization 
gain exhibits a significant coupling when a stabilator failure occurs. However, the 
MIL-F-8785C specification on sideslip to bank angle coupling is that the minimum 
(u n is increased if u> 2 |^/,0| > 20. For our aircraft with a stabilator failure we have 
ul\(f)/P\ < (2.83) 2 (2.2) = 17.6 < 20. Therefore, our responses meet the MIL spec- 
ification even though a smaller sideslip to bank angle coupling may be desirable. 
The responses for the unfailed aircraft to an initial sideslip are shown in Figure 8 
for both the eigenstructure assignment and multiobjective optimization gains. We 
observe that that the sideslip and yaw rate responses are identical and there is only 
an insignificant degradation in the coupling to roll rate. Figure 9 shows the responses 
using both gains for the unfailed aircraft until time t — 0.5 sec when a rudder failure 
occurs. We observe that the eigenstructure assignment gain allows oscilllation and 
large settling time whereas the multiobjective optimization gain yields an acceptable 
response after failure. 
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Figure 5: Yaw rate responses 
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Figure 6: Roll rate responses 


25 




bank angle (deg) 
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eigenstructure assignment gain 
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Figure 7: Bank angle responses 
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